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We study the low energy physics of a Kondo chain where electrons from a one-dimensional band 
interact with magnetic moments via an anisotropic exchange interaction. It is demonstrated that 
the anisotropy gives rise to two different phases which are separated by a quantum phase transition. 

In the phase with easy plane anisotropy, Z 2 symmetry between sectors with different helicity of the 
electrons is broken. As a result, localization effects are suppressed and the dc transport acquires 
(partial) symmetry protection. This effect is similar to the protection of the edge transport in 
time-reversal invariant topological insulators. The phase with easy axis anisotropy corresponds to 
the Tomonaga-Luttinger liquid with a pronounced spin-charge separation. The slow charge density 
wave modes have no protection against localizatioin. 
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Introduction. One-dimensional systems present an 
ideal platform for formation of charge density waves 
(CDW) P; the transport in clean systems is almost ideal 
[ 2 ]. However, for realistic interactions and at low tem¬ 
peratures, even a weak disorder pins the CDW suppress¬ 
ing the charge transport [3]. The ideal transport can be 
protected by symmetries: a well-known example is the 
edge transport in two-dimensional time-reversal invariant 
topological insulators (TIs)[lHZ]. The topologically non¬ 
trivial state of the bulk and time-reversal symmetry lead 
to a lock-in relation between the chirality and the spin 
of edge modes making them helical [5]. As a result, the 
electron backscattering must be accompanied by a spin- 
flip; hence the edge transport becomes immune to effects 
of potential disorder. Other processes which can sup¬ 
press the ideal transport include scattering by magnetic 
impurities [5] or inelastic processes due to interactions 
[IDHIll- All of them become ineffective at low tempera¬ 
tures. The presence of (almost) ballistic edge transport 
has been confirmed in state-of-the-art experiments |15b 
[18]. Hence it is accepted that the ballistic transport is 
protected by time-reversal symmetry and this protection 
is removed when this symmetry is broken iiiiiinj. 

Helical boundary modes can exist in noninteracting 
systems due to topological nontriviality of the bulk pi] . 
In this Letter, we show that helical modes may emerge 
in interacting systems as a result of spontaneous sym¬ 
metry breaking. As an illustration, we study a model of 
Kondo chain p3H3S] consisting of band one-dimensional 
electrons interacting with local spins; the Hamiltonian of 
this system is: 

H = -t'^ -b H.c. (1) 

n m 

Here = (c^(n),c^(n)) are electron operators at lat¬ 
tice site n; (t“ are Pauli matrices (a = x,y,z)-, S°‘(m) 
are components of the spin-s operator located on lattice 
site m; t denotes the overlap integral. It is assumed 


that sites {m} constitute some (not necessarily regular) 
subset of sites {n}. We concentrate on the regime of 
sufficiently high density of spins where the Kondo ef¬ 
fect is suppressed and the physics is determined mostly 
by the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter¬ 
action m- The band is far from half filling, the spins 
are quantum and the coupling constants are much smaller 
than the bandwidth, sJa ^ t. We will consider the cou¬ 
pling which is isotropic in the AY-plane: = Jy = J±. 

Brief summary of the results: The low energy (LE) 
behavior of model 0 includes two distinct regimes cor¬ 
responding to the easy axis (EA), > J_l, and the 

easy plane (EP), < Jj_, anisotropy. In the first 

case, all quasiparticle (fermionic) excitations are gapped. 
The transport is carried by gapless collective modes, the 
charge and the spin density waves. The CDW couples 
to a potential disorder which is able to pin it and to 
block the charge transport. The SU(2) symmetric point, 
Jz — Jx) is the point of quantum phase transition into 
a phase with spontaneously broken helicity. In the EP 
phase at T = 0, quasiparticles with a given helicity ac¬ 
quire a gap and the other helical branch remains gapless. 
The charge transport is carried by the gapless helical elec¬ 
trons and by the slow collective excitations (spin-fermion 
waves). If the spin U(I) symmetry is respected the long 
range helical ordering makes single-particle backscatter¬ 
ing of the gapless modes impossible as in the noninteract¬ 
ing TIs. This leads to suppression of localization effects: 
the localization radius becomes parametrically large and 
the dc transport acquires a (partial) symmetry protection 
in finite but long samples. 

Continuum limit: To describe the LE physics we de¬ 
velop a continuum limit theory. This requires to single 
out smooth modes. We linearize the spectrum of elec¬ 
trons and expand operators c in smooth chiral modes: 

cu{n) = (^) + ^ = <0 ; ( 2 ) 

were ^0 is the lattice constant. The Lagrangian density 
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of the band electrons becomes 




(J 0 i)dr — i(/ 0 T^)vFdx 




(3) 


Here r is the imaginary time; the first space in the tensor 
product is the spin one, the Pauli matrices t“ act in 
the chiral space; / = diag(l,l); vp — 2t^Qsm{kF^o) is 
the Fermi velocity {kp is the Fermi momentum); 4'"'" = 
{R^, is the 4-component fermionic spinor field. 

Contrary to Ref. [22] , where the effects of forward scat¬ 
tering at ~ t (i.e., of the Kondo physics) were consid¬ 
ered, we suggest that the LE physics in the dense limit 
with Ja <C t (dominated by the RKKY interaction) is 
governed by backscattering of the fermionic modes. It is 
described by 

Cks = Ps Y. + (4) 

a—x,y,z m 


Ps denotes the dimensionless spin density. £bs is expected 
to lead to opening of the spectral gaps thus reducing the 
energy of the electrons. As will be clear from the subse¬ 
quent discussion, the resulting physics is quite different 
from that of Ref. [22] . 

We can eliminate the oscillatory factors in @ by ab¬ 
sorbing them into the spin configurations which amounts 
to separation of fast and slow spin variables [55]. The 
standard parametrization of the spin by azimuthal and 
polar angles, S = s{sin(0) cos(i/'), sin(d) sin(^/^), cos(d)}, 
with the integration measure = sm{6)'D{9}'D{ilj} 

[55] is not convenient for our purposes. Therefore, we 
change to the rotating orthonormal basis 61 ^ 2,3 with 
03 = S/s. We define a “longitudinal”, 5||, and the 
“transverse”, iS_l, components of the new spin vector 
S = S± -b 5|| (Fig{^: 



s 


03 sin ay; 


s 


= [01 cosa_L -|-02 sina_L] cosay; (5) 


a_L = 2kF^orn + a{x). The orthonormality can be re¬ 
solved by choosing 


01 = {— cos( 0 ) cos{ijj), — cos( 0 ) sin(^/), sin( 0 )}, ( 6 ) 

02 = {sin(V'), - cos(^/), 0 }. (7) 


The integration measure for a, ay will be Vlfla} = 
cos(ay)I?{ay}I?{a}, the total measure reads = 

'D{Vlot}'D{^s}■ This does not result in overcounting the 
degrees of freedom since we will find a scale separation 
with two fast (massive ay, 0 ) and two slow (massless a, ip) 
angles [30] ■ Verihcation of the scale separation and sta¬ 
bility of the chosen spin configuration will confirm self- 
consistency of our approach. 

Inserting the new parametrization in Eq.Q and keep¬ 
ing only the non-oscillatory terms, we find LE Lagrangian 
Ces = £e + £bP + £wz where 


^(sl) _ sps 

^bs - 2 



J± 



sin^ 



-- e-^^cos^ 





FIG. 1.^ Transformation from the frame of the vector S to 
that of S. Angles define the modulus of the transverse 

—b 

component 5x and its rotation around the longitudinal com¬ 
ponent 5||, respectively. 


-I-Jz sin(0)(j^|Le + H.c.; s = scos(ay); (8) 

£wz is the topological Wess-Zumino term [HIES]: 

£wz = 15 ^ 8 ^ 0 "^ sin(ay)[5ra + cos(0)cIrV’]- (9) 

The fermionic gaps become maximal at ay = 0 and 9 — 
0,7r/2,7r. Thus, we expect three extrema of the action 
whose stability depends on the ratio J±/Jz- 

EA anisotropy, > Jp'. The term 0{Jz) dominates 
and opens the gap in all fermionic modes. This can be 
shown straightforwardly after removing the angles a, pj 
from the backscattering term Q by using the Abelian 
bosonization [SS] [34] : we bosonize the fermions and shift 
bosonic phases: 

l>e = - a/2, 0s = 0. - V’/2. (10) 

Here and 0s are the charge and the (dual) spin phases, 
their gradients are coupled to charge- and spin source 
fields, respectively [3S]: £/, = hcd^^c + hed^Qg. Af¬ 
ter shifting the bosonic phases, terms cx hcda:aj2 and 
cx hsdxip/2 arise in the Lagrangian. Finally, we can re¬ 
turn to the fermionic variables: 

£(si) ~ £g -b +y^ £tl(^’, vf) + £wz- (II) 

2^—a,'Ip 

Here £TL(‘f,'c) = [(5^4?)^ + {v dx^Y]/T^v is the La¬ 
grangian of the Tomonaga-Luttinger Liquid (TLL), i.e., 
the anomaly [36] . 

For fixed values of {9, ay}, the fermionic spectrum con¬ 
sists of the four Dirac modes with the masses given by: 

= (sps/2)^ ^yO^~cos^~0”-bJ^~sin^ ± Jx^ • (12) 

Integrating out the gapped fermions, we get the contri¬ 
bution to the ground state energy: 

Y + oiJl, Jz)- ( 13 ) 

x=± 
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If Jz > J_L, -Bgs has minima at 0 = 7r/2,a|| = 0; small 
fluctuations around the minima read: 


« (J| - J\) cos^{0) + { Jl + Jl_) sin^(a||), (14) 

where £ = \a(t/J){ sPs)‘^^q/A'kvf and we do not dis¬ 
tinguish between Jz and J_l in the logarithm. Using 
Eqs.( 9|14' ) and integrating over the Gaussian fluctuations 
of the angles, we can find parameters of £tl(q:) which 
are renormalized due to the coupling of the spin wave to 
the gapped fermions: Avalvp = Ka ^ 1 [S3- The LE 
Lagrangian for the EA anisotropy is [38]: 

£ea = £tl(V') Vf)/4, + £tl(Q!, Va)/Ka + ■ (15) 


£ea corresponds to two U(l)-symmetric TLL models with 
the slow charge, a, and the fast spin, ip, bosonic modes. 

Breaking symmetry: If Jz 3> J_l, then m+ ~ m_, 
all fermionic modes have (almost) the same gap ~ Jz, 
cf. Eq.([^. Mass m_ progressively shrinks towards the 
SU(2) symmetric point of the quantum phase transition 
where m_ = 0 and one subsystem of the helical fermions 
becomes gapless. Our approach looses its validity at 
TO- —>■ 0. We leave a description of the SU(2) symmetric 
point for future studies and consider instead the case of 
the strong EP anisotropy Jz <C Jj.. 

EP anisotropy, Jz ^ J_l: To make the consideration 
transparent, we put Jz —>■ 0 and rewrite Eq.([^ as a sum 
of helical contributions: 


= sp,JfR\ cos^ (0/2) -b H.c. (16) 

£|^f) = -SpsJfRI sin^ (0/2) + H.c. (17) 

If 0 ~ 7r/2, both helical sectors have a gap though 
the coupling constant Jj_ is effectively decreased because 
sin^ (0/2) ~ cos^ (^/2) — 1/2. If 0 ~ 0,7r, only one heli¬ 
cal sector acquires the gap to = to+( J 2 = 0, ay = 0 = 0), 
and J_L is not suppressed because either sin^ (0/2) — 1 
or cos^ (fi*/2) — 1. Since the contribution of the gapped 
fermions to the ground state energy is negative and 


quadratic in the gap, Eq.(13), we conclude that 0 = 7r/2 
yields maximum of the energy and two (degenerate) min¬ 
ima are 9 — 0,tt. Thus, the Z 2 symmetry between the he¬ 
lical subsystems is spontaneously broken confirming that 
the SU(2) symmetric point is the point of a quantum 
phase transition [39j . 

Let us consider the configuration 0 ~ 0 where only 
yields the femionic gap [30]. One can straight¬ 
forwardly estimate that contributions of the gapped 
and the gapless fermions to fluctuations of the ground 
state energy are of order ^ ( J^/uf) sin^(0/2) and ^ 
{J]_/vf) sin^(0/2), respectively. The latter is subleading, 
it is beyond our accuracy and must be neglected. Thus, 
£(,g is irrelevant for the effective LE theory and must 
be neglected too. The combination ip — a becomes re¬ 
dundant and Ip in the combination ip + a [see Eqs.( 9[T6 )] 
can be absorbed in a: ip + a ^ a m- Now, we can 


proceed very similar to the case of the EA anisotropy: 

(a) eliminate the shifted spin phase a from £/g with 
the help of the transformation 

l>e = + a/2, 0, = 0, - a/2; (18) 


(b) integrate out massive helical fermions and obtain the 
fermionic energy close to its minima: 

5E,^/£ ~ Ji[sin2(0/2) +sin"(a||)/2]; 


(c) integrate out small quadratic fluctuations of angles 
around the stationary value; and (d) bosonize fermions 
from the gapless helical sector by using the Abelian phase 
$H- These steps yield the effective Lagrangian for the 
case of the EP [38] : 

£ep = £TL(‘i>H,t'F)/2 + £TL(a,<)/i^; + £|r'’^(19) 


where Av'^/vf = K'^ ^ 1 [37] • Similar to the EA 
anisotropy, £ep corresponds to two U(l)-symmetric TLL 
models with the fast, <1 >h, and the slow, a, bosonic modes. 
However, as we discuss below, the effective theories with- 
and without the helical symmetry have different trans¬ 
port properties if a disorder is added. 

To conclude this section, we note that Eqs.( IGjlg ) are 
equivalent to their counterparts describing a helical edge 
mode in the TI with an array of the Kondo impurities [191 
[3n| . In our case, however, this helical mode has emerged 
as a result of spontaneous symmetry breaking. 

Density correlation functions and effects of the dis¬ 
order. Let us keep source terms for the charge sector: 

= /ic9xa/2; £j/P^ = /ic(9a;4>H+9a;a/2). The charge 
density-density correlation function at given frequency 
and momentum reads as: 

Cea oc Q"(a*a); Cep oc + (a*a)/4) . (20) 


Cea,ep with Lagrangians £ea,ep correspond to the ideal 
metallic transport. In the EA case, it is supported by 
the slow CDW with the small compressibility, ATq,. Cep 
contains the contribution from the helical quasiparticles 
with the bare velocity and from the slow collective wave 
with the small compressibility, K'^ [33]. 

The coupling of backscattering spinless impurities to 
the fermions is described by: 


Vdis[g] = H.c. (21) 


Here g{x) is the smooth 2fci?-component of the scalar ran¬ 
dom potential. We use the model of the Gaussian white 
noise: (5^’^)dis = 0; ( 5 (a:i) 5 *(a; 2 ))dis = D6{xi - X 2 ), as¬ 
suming that the disorder is weak, D <C (to±, m)vF, and it 
cannot change the gaps. After shifts Eq.( 10|18), the po¬ 


tential g acquires the phase factor: g ^ g x Thus, 

the backscatterering impurities are coupled to all gapless 
charge carriers (collective waves and helical fermions). 

To figure out whether the disorder may lead to local¬ 
ization, we perform the disorder averaging and integrate 
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out the massive fermions [43) . The relevant terms ap¬ 
pear only in P^-order and have a different form in EA 
and EP phases. In the Hrst case, V'^ couples directly to 
exp(iQ;); in the EP phase, it couples to i?+L_cr exp(ia). 
The latter fact is related to impossibility of single par¬ 
ticle backscattering in the phase with broken helicity. 
The power counting indicates the parametric difference 
in the localization radius in both phases: ~ 

Kj^Djvp'rnf'l^ <C 1, with ~ {vFlrn) {vpinn/'D)^. 


Localization can block the dc transport if a sample 
size L is large: L ^ We thus conclude that the 

ballistic transport in the phase with broken helical sym¬ 
metry acquires the symmetry protection up to the para¬ 
metrically large scale This conclusion holds true 

as long as the U(l) symmetry in the spin sector is re¬ 
spected. Breaking the U(l) spin symmetry (e.g. after 
introducing an anisotropy in the XY-plane) allows the 
direct backscattering of all fermions and removes protec¬ 
tion of the ideal transport in the EP phase [cf. local¬ 
ization of the helical edge modes of the TIs m in the 
absence of the U(l) spin symmetry]. 


Finite temperature effects in the clean case: All pre¬ 
vious calculations have been done for zero temperature, 
T = 0. They can be generalized for T ^0 provided 
T is smaller than the fermionic gaps. Finite temper¬ 
ature restores a broken helical symmetry at Jz < J_l 
since thermal fluctuations produce domains with oppo¬ 
site helicity. When the spin configuration interpolates 
between the phases with different helicity there is an en¬ 
ergy increase of the order of the difference between the 
energy in the unstable state (with 0 ~ 7r/2) and the 
energy of one of the ground states (with 0 ~ 0 or tt). 
Thus, we can estimate the energy of the domain wall as 
Awaii cf.Eq.(13). The maximal number of 

the domain walls in the system of the size L reads as 
Lm/vF- If T ^ Ewaii, it becomes exponentially sup¬ 
pressed: A^wall ^ Lm/n_Fexp(-Awall/T)- If A^wall > 1, 
the walls appear and block the quasiparticle transport 
since the electrons with a given helicity are massless 
only in one domain and massive in the other (neighbor¬ 
ing) one. Hence the electrons are reflected from domain 
boundaries. On the other hand, an influence of the do¬ 
main wall on the field a is reduced to a jump in the 
Luttinger parameter K'^ which cannot affect the dc con¬ 
ductance, cf. Ref. [13]. Thus, we arrive at the conclu¬ 
sion that the dc transport in the phase with the broken 
helical symmetry will remain ballistic even at finite tem¬ 
peratures. Temperature effects in the disordered case de¬ 
serve a separate study because of a complicated interplay 
between formation of the domain walls and many-body 
(de)locahzation of collective waves [15] . 

Validity: The effective LE theory, Eqs.( 15|l9 ), is valid 
at energies below the smallest fermionic gap, to_ and 
m for the EA and the EP anisotropy, respectively. Since 
TO- vanishes at the SU(2) symmetric point, the approach 


fails in the vicinity of the quantum critical point. Quickly 
oscillating contributions cx which we neglected, 

are generically unable to change the physics at the large 
distances: If the Kondo chain is close to incommensura¬ 
bility the quickly-oscillating exponentials can be treated 
as random variables, cf. Ref. [53]. We note, however, 
that, in the most interesting case of the broken helical 
symmetry, the amplitude of the oscillating terms is sup¬ 
pressed in the vicinity of the classical spin configuration, 
0 ~ 0, as ~ (^O'^I/'^^f) sin'^(0/2) [see the discussion of 
the derivation of Eq. @] which is squared after averag¬ 
ing over the random fluctuations, i.e., becomes negligible. 

Conclusions: We have demonstrated that the dc charge 
transport in the Kondo chain model Q with the U(l) 
symmetry of spins remains ballistic in long samples, 
L < even in the presence of the potential disor¬ 

der when the anisotropy of the exchange interaction is of 
the easy plane type. Due to the spontaneous breaking 
of the Z 2 symmetry the current is carried by quasiparti¬ 
cles possessing a particular helicity (i.e. whose spin and 
chirality are locked) and by composite spin-fermion col¬ 
lective modes. In the presence of the U(I) spin symmetry, 
all gapless modes are protected from simple backscatter¬ 
ing by the mechanism similar to that in noninteracting 
TIs. We emphasize that the symmetry protected trans¬ 
port in our model results from interaction many-body 
effects instead of the coupling to the non-interacting and 
topologically non-trivial bulk. In the case of the easy 
axis anisotropy, the helical symmetry is respected. The 
quasiparticles are fully gapped and the transport is car¬ 
ried solely by the collective modes. The slow CDWs do 
not posses the symmetry protection: the potential disor¬ 
der can pin them and render the Kondo chain insulating. 

A.M.T. acknowledges the hospitality of Ludwig Max¬ 
imilians University where this work was done. A.M.T. 
was supported by the U.S. Department of Energy (DOE), 
Division of Materials Science, under Contract No. DE- 
AC02-98CH10886. O.M.Ye. acknowledges support from 
the DEG through SFB TR-12, and the Cluster of Excel¬ 
lence, Nanosystems Initiative Munich. We are grateful 
to Vladimir Yudson, Igor Yurkevich for useful discussions 
and to Dennis Schimmel for carefully reading the paper 
and for his participation in the derivation of the Wess- 
Zumino term. 


[1] T. Giamarchi, Quantum physics in one dimension 
(Clarendon; Oxford University Press, Oxford, 2004). 

[2] A. Rosch and N. Andrei, Phys. Rev. Lett. 85, 1092 

( 2000 ). 

[3] T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 
(1988). 

[4] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

( 2010 ). 

[5] X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 





5 


( 2011 ). 

[6] S.-Q. Shen, Topological insulators: Dirac Equation in 
Condensed Matters (Springer, 2012). 

[7] M. Franz and L. Molenkamp, Topological Insulators (El¬ 
sevier Science, 2013). 

[8] C. Wu, B. A. Bernevig, and S. C. Zhang, Phys. Rev. 
Lett. 96 , 106401 (2006). 

[9] Y. Tanaka, A. Furusaki, and K. A. Matveev, Phys. Rev. 
Lett. 106 , 236402 (2011). 

[10] T. L. Schmidt, S. Rachel, F. von Oppen, and L. 1. Glaz- 
man, Phys. Rev. Lett. 108 , 156402 (2012). 

[11] V. Cheianov and L. 1. Glazman, Phys. Rev. Lett. 110 , 
206803 (2013). 

[12] J. 1. Vayrynen, M. Goldstein, and L. 1. Glazman, Phys. 
Rev. Lett. 110 , 216402 (2013). 

[13] J. 1. Vayrynen, M. Goldstein, Y. Gefen, and L. 1. Glaz¬ 
man, PRB 90, 115309 (2014). 

[14] N. Kainaris, 1. V. Gornyi, S. T. Carr, and A. D. Mirlin, 
Phys. Rev. B 90 , 075118 (2014). 

[15] M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X. L. Qi, and S. C. Zhang, 
Science 318, 766 (2007). 

[16] A. Roth, C. Briine, H. Buhmann, L. W. Molenkamp, 
J. Maciejko, X.-L. Qi, and S.-C. Zhang, Science 325, 
294 (2009). 

[17] 1. Knez, R.-R. Du, and G. Sullivan, Phys. Rev. Lett. 
107 , 136603 (2011). 

[18] K. Suzuki, Y. Harada, K. Onomitsu, and K. Muraki, 
Phys. Rev. B 87, 235311 (2013). 

[19] B. L. Altshuler, 1. L. Aleiner, and V. 1. Yudson, Phys. 
Rev. Lett. Ill, 086401 (2013). 

[20] O. M. Yevtushenko, A. Wugalter, V. 1. Yudson, and 
B. L. Altshuler, “Transport in helical Luttinger Liquid 
with Kondo impurities,” (2015), arXiv: 1503.03348. 

[21] B. A. Bernevig and T. L. Hughes, Topological insula¬ 
tors and topological superconductors (Princeton Univer¬ 
sity Press, 2013). 

[22] O. Zachar, S. A. Kivelson, and V. J. Emery, Phys. Rev. 
Lett. 77 , 1342 (1996). 

[23] G. Honner and M. Gulacsi, Phys. Rev. Lett. 78 , 2180 
(1997). 

[24] H. Tsunetsugu, M. Sigrist, and K. Ueda, Rev. Mod. 
Phys. 69 , 809 (1997). 

[25] E. Novais, E. Miranda, A. H. Castro Neto, and G. G. 
Cabrera, Phys. Rev. B 66, 174409 (2002). 

[26] M. Gulacsi, Adv. Physics 53, 769 (2004). 

[27] This statement means that the spin density is sufficiently 

high so that the Kondo effect is cut by the gap generated 
by backscattering terms in H, cf. Refs. [T^ . ft holds 

true when the gap exceeds the Kondo temperature. 

[28] This approach is supported by numerical studies of the 
Kondo chain of classical spins which indicate that the 
phase diagram of the system reflects helical spin configu¬ 
rations [4^ . The physics, which is governed by the helical 
effects in spins, is also studied in R.efs. [47(449) . 

[29] The topological Wess-Zumino term (the Berry phase), 
£wz, must be added to the Lagrangian m- 

[30] We note that the procedure of separation of slow and 
fast variables is standard and is routinely used in RG 
calculations: see, for example, Chapter 1 in Ref. m- 

[31] A. M. Tsvelik, Quantum Field Theory in Condensed Mat¬ 
ter Physics (Cambridge: Cambridge University Press, 
2003). 

[32] Derivation of Eq. is explained in Suppl.Mat. No.l. 


[33] 

[34] 

[35] 

[36] 


[37] 

[38] 


[39] 


[40] 

[41] 

[42] 

[43] 

[44] 

[45] 

[46] 


[47] 

[48] 

[49] 

[50] 


A. O. Gogolin, A. A. Nersesyan, and A. M. Tsve¬ 
lik, Bosonization and strongly correlated systems (Cam¬ 
bridge: Cambridge University Press, 1998). 
Alternatively, one can do gauge transformations of the 
fermionic fields and obtaine an anomaly from the Jaco¬ 
bian 

Two remaining fields ca n be integrated out. 

While deriving Eq.(lll, we have neglected the gra¬ 
dient coupling of the fermions and the spin phases, 
T(a,'!/')]4', which is justified for energies below the 
fermion gaps. 

Expressions for the Luttinger parameters and K'^ are 
given Suppl.Mat. No.3. 

While deriving Eqs.(15l and (191, we have neglected 
in /Iwz the subleading terms oc sin(Q;||) cos(@)9rV' and 
oc sin(a||) sin^(0)9T'!(’, respectively. This is justified for 
energies below the fermion gaps. 

We note that helicity is given by the product of the 
spin projection and the electron velocity, H = sign(cr?;). 
Hence, helical symmetry is discrete and it can be spon¬ 
taneously broken in one dimension at T = 0. The cor¬ 
responding order parameter is related to the average 
Ca = {Sb{x)Sc{x-\-^o)^abc)- iu the phase with the broken 
symmetry Cz oc sin(2fc_F^o) aud Cz = 0 otherwise. 

The second minima 6 ~ tt can be considered simi larly 
and it yields the same effective Lagrangian £ep, Eq.(19l. 
One can arrive at the same conclusion even faster using 
the standard parametrization of S, see Suppl.Mat. No.2. 
The small compressibility of the bosonic modes leads to 
suppression of the Drude weight which reflects the cou¬ 
pling of the spin waves to the gapped fermions [Si- 
All essential details of this calculus are presented in 
Suppl.Mat. No.4. 

N. Sedlmayr, J. Ohst, 1. Affleck, J. Sirker, and S. Eggert, 
Phys. Rev. B 86, 121302 (2012). 

D. M. Basko, 1. L. Aleiner, and B. L. Altshuler, Annals 
of Physics 321, 1126 (2006). 

W. Hu, R. T. Scalettar, and R. R. P. Singh, “In¬ 
terplay of magnetic order, pairing and phase separa¬ 
tion in a one dimensional spin fermion model,” (2015), 
arXiv: 1506.04809. 

B. Braunecker, P. Simon, and D. Loss, Phys. Rev. B 80, 


165119 (2009). 

J. Klinovaja, P. Stano, A. Yazdani, and D. Loss, Phys. 
Rev. Lett. Ill, 186805 (2013). 

T.-P. Choy, J. M. Edge, A. R. Akhmerov, and C. W. J. 
Beenakker, Phys. Rev. B 84, 195442 (2011). 

A. Grishin, I. V. Yurkevich, and I. V. Lerner, Phys. Rev. 
B 69, 165108 (2004). 






6 


SUPPLEMENTAL MATERIALS 
1. Derivation of the Wess-Zumino term 


Here we discuss the subject of spin action which is usually formulated in the Wess-Zumino form m- This form is 
invariant under rotations, however, it requires an integration over an auxiliary variable which is not convenient for our 
purposes. We have to find another formulation which would not include the additional integration and, nevertheless, 
would allow us to change the basis. 

Let us start with the spin being defined as S = s 63 , see the paragraph before Eq.([^ in the main text. Here 63 is one 
of vectors from the orthonormal basis 01 ^ 2 , 3 , an example of the vectors 61^2 is given in Eqs.(6|7). The Wess-Zumino 
term for this representation is well-known: 


>Cwz[6',V'] = cos(6>)5^V'• 


( 22 ) 


The boundary contribution has been neglected in Eq.(22) and in all equations below since we are interested in smooth 
(semiclassical) spin modes. We note that Eq.( 22) is invariant with respect to all 0(2) rotations of the vectors ei 2 T 03 : 

ei = cos(,5)ei -|- sin(; 3 )e 2 , ei = - sin(/3)ei -|- cos(/ 3 )e 2 . (23) 


This is because (3 is the gauge angle and it either drops out from Eq.(22| (if /3 = const) or yields an unimportant 
boundary contribution (if depends on time). This allows us to rewrite Eq.(22) by using the vectors 63^2 from any 
orthonormal basis: 


C.wz[OA] = ^[(e 2 ,i 9 ^ei) - (01,9^62)] = ^(e 2 ,i 9 rei). 


(24) 


The first equality in Eq.(24) can be verified by direct inspection after inserting expressions ( |6|7[ ) into Eq.([^ and the 
second one follows from the orthogonality condition ( 01 , 62 ) = 0. We note that Eq.(|24|) contains only scalar products 
of two vectors and, therefore, it is invariant under the global rotation of the {a;, y, z}-basis: 


— (Ej, C^t-E/j;) if El 2 — R-xyz ^ 1 ^ 2 ^ 

To avoid confusions, subscriptips of orthogonal matrices show the basis where they operate. 
Now we change to the new spin, Eq.®, and define the new orthonormal basis 03 23 with: 


(25) 


63 = 5/s = cos(Q!||)[cos(a_L) 6 i -I- sin(Q!_L) 62 ] -I- sin(a ||)63 . 


(26) 


Two remaining vectors from the new basis can be chosen, for example, as follows: 

(27) 

(28) 


0 / = — sin(a||)[cos(a_L) 6 i -I-sin(a_L) 62 ] -I-cos(a||) 63 , 
62 = sin(aj_) 6 i — cos(q:j_) 62 . 


Let us first assume that aj_j| does not depend on time. In this simple case, the transformation (26-281 is global in 
the 63 2,3 basis but it is local in the {x, y, zj-basis. The latter statement results from the rotation oFtne basis vectors 
Gi, 2 , 3 - Thus, before using Eq.(26-28), we have to rewrite Eq.(24) in the form which is invariant under all possible 
global rotations. Such a form reads as 


(S, 6i) (Oj-, 


(29) 


where Cijk is the antisymmetric tensor. Eq.(29) reduces to Eq.(24) if S = 503 . It is clearly invariant under the global 
rotation by the matrix i?xyz- The invariance under the global rotation by the matrix i?i 23 can be easily shown after 
rewriting the vector S in the 63 ^ 2 , 3 -basis, substituting the basis (ej, 63 , 03 )"^ = 7?i23(ei, 62 , 63 )"'" into Eq.(29) and using 
the identity 


E 


[-Sl23]ip [-^123]jm ['^123]A;g ^ijk 


2j,/c=1,2,3 


^pmq • 
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Following the approach explained in Sect. Continuum limit of the main text, we insert Eqs.(26-28) into Eq.(29) 
allowing the angles a_L,|| to be the independent variables and substitute s sin(a||)e 3 for the spin S. This procedure 
yields Eq.([^ in the main text. The second substitution takes into account an effective decrease of the slow component 
of the spin: after introducing the new rotating frame, the size of the spin will be given by the overlap of this effective 
spin with the original one, i.e. by 5 ( 03 , 63 ) = ssin(Q:||), see Fig 0 

To finalize the discussion of the Wess-Zumino term, we point out a short-cut which allows one to obtain the answer 
Eq.® even faster: we can 1) exploit the gauge invariance described in Eq.(23) and directly insert the vectors e'^ 2 
into^q.(24); 2) allow the angles a_Lj| to be independent variables in the actions; 3) do the shift of a± and omit the 
oscillatory terms in £wz- 


2. Alternative derivation of the effective Lagrangian in the case Jz = 0. 


Let us consider the extreme case of the easy plane anisotropy where = 0. Eq. simplifies to 

/lbs = p. {s+rIl^ + S-R\l^') + H.c. (30) 

here = 5'“ ± . Now we can use the standard parametrization of S by azimuthal and polar angles: 

S = s[sin(P) cos(A), sin(P) sin(A), cos(P)], (31) 

^ Pbs = spsJ± cos(P)e2i'=^^ -f H.c. (32) 

(notation for the angles are changed as compared to the main text to avoid confusions) and note that the slow spin 
modes can be easily singled out after a shift 


A^A±2kFX. (33) 

The choice of the ± sign substitutes now the choice of the minima (either 0 = 0 or 0 = tt, see the main text): it 
breaks the helicity symmetry and corresponds to the strong correlation of the spins to one or the other sector of the 
helical fermions. The helical sector, which is not correlated with the spins, acquires dfci^’-oscillations and vanishes in 
the effective Lagrangian for the backscattering. For example, choosing the plus sign, we obtain 

) = sp^Jf cos(P)e-'^P|L; + H.c., ^ = 0; (34) 

with the Wess-Zumino term £wz = isPsCo"^ cos{P)drA and with the integration measure P{n} = sin(P) 'D{P}'D{A}. 
This confirms that the third angle becomes redundant for the effective low-energy theory at Jz = 0. 


3. Luttinger parameters Ka and K'^ 


Calculations described after Eq.(14| of the main text yield the following expressions for the Luttinger parameters: 

- £n \/log(t/J) < 1; K'^ \/log(t/w) < 1. (35) 


ttvf 


TTSPsVF 


4. Derivation of the localization radius, 

In this section, we derive an estimate for the localization radius in the Kondo chain coupled to spinless backscattering 
impurities. Firstly, we replicate fields 

a G 

and calculate the Gaussian integral over the random field g. This yields the standard contribution to the action 

{R^L)'^ [1] X (e-i“/"(P^P)) [2] + H.c ., n={x, r„, a„, a„} , n = 1, 2 . (37) 

0 ^ 1,2 <^1,2 
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At the next step, we integrate out massive fermions perturbatively by doing an expansion in the small parameter 
'D/{m±^m)vF <C 1. Our goal is to find leading terms which can result in pinning of all collective charge carriers 
(EA and EP phases) and in localization of the massless helical fermions (EP phase). We do it separately for the two 
different phases. 


4-1 The EA phase: most relevant terms 


Let us put Jj_ —)■ 0 and introduce a spin-dependent fermionic mass mea,{<j) = ±(sps/2) Jz where the plus (minus) 
correspond to a =t (cr =),). This allows us to simplify the derivation without loss of generality. The matrix Green’s 
function for the fermions with a given spin reads: 





^ -mea{cr) 


-mea(o-) 




(38) 


where G^^\ are the Green’s functions of free chiral particles. It is important that Gm is short ranged and it decays 
beyond the time scale Ijmes. (or beyond the coherence length ,^ea = Vf/rriea.)- 

Leading terms of the order of 0(T>^) are given by (S'dis) where brackets mean that the massive fermions are 
integrated out. The corresponding diagrams are shown in Fig[^ It is easy to check that the diagrams from Figj^a 
cancel out after summation over spins indices because mea(t) = ~Wea(4-). The diagrams from Figj^b are trivial since 
Gm is diagonal in the replica space and the spin phase a is smooth on the scale l/wea; therefore, 


gia[l]/2g-ia[2]/2 _ ^ ^ 


(39) 


with some small gradient corrections which are unable to yield pinning. 



FIG. 2. First order diagrams 0{V^) for the EA phase. Red (green) triangulars denote e'“f^ (e“*“^^) with argnments of 
either the 1st or the 2nd vertex; dashed lines are the disorder correlation fnnctions, solid lines stand for Green’s functions of 
the massive fermions. 

Sub-leading terms of the order of 0{'D‘^) are given by (5'dis'S'dis) which generates a lot of diagrams. We leave a 
detailed analysis for the future and calculate only one typical diagrams which survives after all summations and is 
able to generate pinning. An example of such a diagram is shown in Figj^ 



FIG. 3. A typical non-trivial diagram, of the order 0{'D^) for the EA phase; notations are explained in the caption of 

Fig|l 


Neglecting unimportant numerical factors, the analytical expression for the diagram from Figj^ reads as: 


/ d{a;,a;';'ri.2,T( 2}e 


l(a[l]-a[2]) 


G„(l,l') G^(l',l) G^(2,2') Gmi2',2) 


1:2 


(40) 
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Here, we have taken into account the the diagonal structure of Gm results in oi = a'i,a 2 = a!^ and fused together 
slow spin phases, for instance: a[l] + Q![l'] — 2a[l]. Now we note that Gmi^, 1^ = Gm(l — 1^ ^-nd integrate over all 
primed variables: 




i 




f ^ 


\VFmea 


)■ 


(41) 


The structure of Eq.(41) corresponds to the disordered Sine-Gordon model which appears in the theory of the usual 

(42) 


TLL [T] . The effective disorder strength V is renormalized and obeys the well-known RG equation [3] : 

EH: 9ioglog(P)=3-2K„~3, i>{U) = i’o; 


the second equality of Eq.(42) has been obtained by using Eq.(35|. 


4-2 The EP phase: most relevant terms 


We start again from the leading diagrams generated by (S'dis)- Tbe principal difference of the EP phase from the EA 
one is that the matrix Green’s function, Eq.(38), corresponds now to the massive fermions with a given helicity. This 
changes the structure of the first order diagram, see Eig|^ All these diagrams correspond to forward-scattering of the 
massless helical fermions and they contain only small gradients of the phase a, cf. Eq.(39) and its explanation. Thus, 
the leading diagrams are trivial and they cannot yield localization, the sub-leading diagrams must be considered. 



12 12 


FIG. 4. Two typical examples of first order diagrams 0('D^) for the EP phase. Incoming red (green) arrows denote the 
product of smooth fields Ri) with arguments of either the 1st or the 2nd vertex; outgoing arrows denote the 

conjugated product, dashed lines are the disorder correlation functions, solid lines stand for Green’s functions of the massive 
helical fermions. 


(iS'disS'dis) generates 16 diagrams with back-scattering of the massless fermions and exponentials of the phase a 
which do not cancel, see FigNeglecting unimportant numerical factors, the analytical expression for the diagram 


FIG. 5. 
Fig@ 



A typical non-trivial diagram, Dep\ 


of the order 0(T>^) for the EP phase; notations are explained in the caption of 


from Figj^ reads as: 


Gm(2,2') 


(43) 


see explanations after Eq.(40) and note the m must be substituted for mea(o’) in Gm- Galculating integrals over all 
primed variables, we find: 

/ V 


DgJ^xVoY / cl{x;ri,2}e'(“W-«[2])it[2]i?t[i]^ji]^^[2], Vo = V 


\VFm J 


(44) 
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This equation also can be reduced to the form of Eq. (41) if remaining fermions are bosonized and we explicitly single 
out new charge- and spin- density waves. However, the RG equation for T) can be obtained without such a complicated 
procedure with the help of the power counting. Firstly we note that the scaling dimension of each back-scattering 
term in Eq.(44), R and R!^L, is 1. The anomalous dimension of each exponential, e^^“, is K'^ 1. The normal 

dimension in Eq.(44| is 3 which comes from three-fold integral. Combining these dimensions together and neglecting 
small we find 


EP : Slog log(P) =3-2x1 + O(K^) ~ 1; = Vo , ^ep = vp/m . (45) 


4-3 Comparison of the localization radius in different phases 


The solution of the RG equations, Eqs.(42|4'5 ), reads as 


V{x) = Vq 


, V{x) = Vq 




(46) 


with ^ep = vpjm. The localization radius is defined as a scale on which the renormalized disorder becomes of the 
order of the cut-off: 


V 


= K^vl/U Klvl/i,. 


V 


= 4/?ep . 


(47) 


The additional small factor is the equation for Leacan be justified with the help of the standard optimization 
procedure [T] where is defined as a spatial scale on which the typical potential energy of the disorder becomes 

equal to the energy governed by the term oc {dxa)^ in the Lagrangian Lea, Eq.(15|. 

Definitions Eq.(47l result in 


^ea Sea.f'.Q 


^eaVo 


1/3 


UKc 


/nFmea\2/3 


V V 


-) 


r (loc) 
ep 


Vo 


V ) 


(48) 


Assuming ^ea ~ ■fep and Wea m, we obtain 


L 

L 


(loc) 

ea 


(loc) 

ep 


Kr 


(—) 

\Vpm J 


4/3 


< 1 . 


(49) 


This demonstrates that the strong suppression of localization can occur in the EP phase where the helical symmetry 
is broken. 

We note in passing that the scaling exponent of V{x) is the same as in the case of non-interacting Id fermions but 
suppression of localization in the EP phase is reflected by the additional large factor vpm/V in the expression for the 
localization radius L^p. 













